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ABSTRACT 


The calculation of 550 km X 550 km mean gravity anomalies from 1° X 1“ mean 
free-air gravimetry data is discussed. The block estimate procedure developed by 
Kaula is used. Estimates for 1452 of the 1654 blocks are obtained. 


RESUME 


Dans cet article on discute du calcul des anomalies de la gravite 
moyenne sur 550 km x 550 km a partir de donn^es de gravim€trie moyenne 
sans air sur 1° x 1°. On utilise le procede d’ estimation de blocs deve' 
loppie par Kaula. On obtien des estimations pour 1452 des 165 1 * blocs. 


KOHCIIEKT 

OOcyjfCflaeTCH bliumc neHne cpeaHitx rpaBHTau wohhhx aHOMa.JiMw 
550km X 550km, kcxobh m3 cpeflHMX naHHtix rpaBHTOMeTpHH cbo<5ohhcm 
aTMOcfiepbi 1° X 1°. XnoTpetfnHeTCH npouenypa oueHKB OnoKa pa3pa- 
(SoTaHHaa Kayna. IIojiyueHbi oqeHKH hjih 14-52 m3 1654- Sjtokob. 
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THE ESTIMATION OF 550 km X 550 km 
MEAN GRAVITY ANOMALIES 


M. R. Williamson and E. M. Gaposchkin 


1. THE BLOCK COVARIANCE METHOD 


The objective is to obtain estimates of mean gravity anomalies for 550 km X 550 km 
regions from an incomplete set of 1° X 1° mean free-air gravity anomalies. The esti- 
mate procedure is based on the covariance analysis developed by Wiener (1966), 
Kolmogoroff, and Kaula (1967). A description and evaluation of the procedure are 
given in Gaposchkin (1973). 


To obtain the mean gravity anomaly for a 550 km X 550 km block, estimates of all 
110 km X 110 km unit mean gravity anomalies in the block are needed. For units 
where measurements of 1° X 1° mean gravity anomalies are available, the unit mean 
anomalies are taken to be the average of the 1° X 1° mean anomalies within the unit. 
The other unit mean anomalies g are calculated from the measured unit anomalies f 
in the same block by the equation 



t 


where K., 
Jk 


are elements of the block covariance matrix, which is given by 


K jk =K(f ’V • 


K is the intra-block covariance function of the unit anomalies f, and t.^ is the distance 
between the jth and kth units. The covariance function is estimated from 


This work was supported in part by Grant NGR 09-015-002 from the National 
Aeronautics and Space Administration. 
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The block anomalies are obtained by averaging the unit anomalies. The mean anomalies 
for blocks that include no measured 1° X 1° mean anomalies cannot be estimated by this 
procedure. 
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2. THE 1° X 1° DATA 


We obtained four sets of 1° X 1° mean free-air gravity anomalies. A set of 29, 209 
measured 1° X 1° means was obtained from the Defense Mapping Agency Aerospace 
Center (DMAAC) (1973). A set of 1454 1° X 1° means for Australia was obtained from 
Mather (1970). Two sets of 1° X 1° means were obtained from Talwani: 4250 1° X 1° 
means for North America and the North Atlantic (Talwani, Poppe, and Rabinowitz, 
1972), and 3944 1° X 1° means for the Indian Ocean (Kahle and Talwani, 1973). The 
sets of data were combined, the Mather and the Talwani data were used for regions 
when available. The combined set has measured values for 31, 654 of the 64, 800 1° X 1° 
areas. Figure 1 is a map showing the distribution of the data. 

The DMAAC data were compared with the other sets of data and with a previous 
compilation by the Aeronautical Chart and Information Center (ACIC) (1971) at common 
points. Table 1 shows the results of these comparisons. 


Table 1. Comparison of 1° X 1° mean gravity anomalies with DMAAC (1973). 


Data 

Number of 
points compared 

Mean difference 
(mgal) 

RMS 

(mgal) 

Australia 
(Mather, 1970) 

1364 

1. 64 

24. 16 

North America and 
the North Atlantic 
(Talwani et al. , 1972) 

3613 

-0. 18 

15.29 

Indian Ocean 
(Kahle and Talwani, 
1973) 

2226 

-1. 66 

23.09 

Worldwide 
(ACIC, 1971) 

19, 164 

-0.23 

16. 99 
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3. THE BLOCK GRAVITY ANOMALY ESTIMATES 


The estimates were obtained with the procedure described in Section 1 by using 
the composite data set described in Section 2. 

The reference ellipsoid used in the calculation is the Geodetic Reference System 
of 1967 (International Union of Geodesy and Geophysics, 1967). A reference ellipsoid 
is determined by the geocentric gravitational constant GM, the rotation rate of the 
earth to, the semimajor axis of the earth a, and the dynamical form factor Jg» The 
1967 system is defined by the ellipsoid 

20 3 -2 

GM = 3. 98603 X 10 cm sec , 

-6 -1 

to = 7. 292 115 1467 X 10 rad sec , 
a = 6. 378160 X 10 8 cm , 

J 2 = 0.001 0827 . 

From these quantities, the flattening f is determined to be 
f= 1/298.247 167 427 . 

A formula for the reference gravity y as a function of latitude <\> is 

o 4 

^1.967 = 978031 ‘ 85 (1 + °* 005 278 895 sin 4> + °* 000 023 462 sin «N m £ al • 

This formula has an accuracy of 0. 004 mgal. Another formula, which has an accuracy 
of 0. 1 mgal, is 

•y = 978031. 846 (1 + 0. 005 302 86 sin 2 4> + 5. 82 X 10 -6 sin 2 2t}>) . 

iyo i 

The ACIC data were compiled using this reference system and the International Gravity 
Standardization Network, 1971 (Morelli and Gantar, 1974). The other data, however, 
are referred to the 1930 international gravity system and the Potsdam Network. To 
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convert these data to the 1967 system, a correction A must be added to the anomalies. 

The 1930 international gravity system is defined by the reference gravity 

^ 1930 = ^78049. 0 {1 + 0, 005 2884 sin^ <j) + 5. 9 X 10 ^ sin^ 2<(>) mgal . 

This reference gravity corresponds to an ellipsoid defined by (Heiskanen and Moritz, 

1967) 

90 9 -9 

GM = 3. 986 329 X 10 cm sec , 

« = 0.729 211 510 X 10 -4 rad sec -1 , 
a = 6.378 388 X 10 8 cm , 
f - 1/297 . 

With the inclusion of the Potsdam correction of 14 mgal, the correction A is 

Y 1930 " Y 1967 “ 14 * 

To an accuracy of 0. 1 mgal 

A= 3. 14 - 13. 58 sin" tj> + 0.02 sin^ 2c|> mgal . 

This formula was used to convert the 1° X 1° mean anomalies to the 1967 system. If 
the anomalies are to be combined with satellite data, the effect of the atmosphere must 
be included. The surface gravity given by the anomalies is increased by 0. 87 mgal to 
obtain the total gravitational force that influences satellite motion. This correction has 
not been made to the anomalies given in this report. 

The covariance function is given in Figure 2 and Table 2. The covariance matrix 
is given in Table 3, The composite set of 31,654 1° X 1° mean anomalies provided 
23, 777 of the 41, 350 unit mean anomalies. From these unit mean anomalies, the esti- 
mate procedure provided 1452 of the 1654 block anomalies. The block anomalies referred 
to the 1967 system are given in Table 4. A punched card deck can be provided on 
request. 
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Figure 2. The block covariance function of unit gravity anomalies. 

We compared these block anomalies with 1183 block anomalies calculated previously 
(Gaposchkin, 1973). The mean difference is 0. 6 mgal and the RMS difference is 
10 mgal. 
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Table 2. The block covariance function of unit gravity anomalies 


Average angular 
distance 

Covariance function 
(mgal 2 ) 

0° 

948 

0.29 

1614 

0.94 

607 

1.21 

457 

1.79 

377 

2. 18 

316 

2. 81 

279 

3. 16 

252 

3.69 

243 

4.20 

229 

4. 77 

182 

5. 22 

161 

5.67 

223 

6.20 

118 

6.69 

201 
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Table 3. The block covariance matrix. 





1 

2 

3 

4 

5 

6 

7 

8 

9 

JO 

11 

12 

13 

1 

9. 4626E+02 

5.717SE+02 

3.4442E+02 

2. G393E+G2 

2. 3447 E+ 02 

5.7175E+02 

4. 2864E+02 

3. 1291E+02 

2. 5177E+02 

2. 3121E+02 

3.4442E+02 

3. 1291E+02 

2.7692E-02 

2 

5. 7175E+02 

9.4826E+02 

5. 7175E+02 

3.4442E+02 

2. 6393E+02 

4, 2863E+02 

5, 7175E+02 

4. 2863E+02 

3. 129 IE+02 

2. 5177E+02 

3. 1291E+02 

3.4442E*02 

3. 1291E+02 

3 

3.4442E+02 

5.7175E+02 

9, 4826E+02 

5. 7175E+02 

3. 4442E+02 

3. 1291E+02 

4.2863E+02 

5.7175E+02 

4, 2863E+02 

3. 1291E+02 

2. 7692E^02 

3. 1291E+02 

3.4442E+02 

4 

2. 6393E+02 

3.4442E+02 

5. 7175E+02 

9. 4 82 6 E+ 02 

5- 7175E+02 

2. 5177 E+ 02 

3. 1291E+02 

4.2863E+02 

5. 7 175E+02 

4. 2863E+02 

2.4411E+02 

2.7692E+02 

3. 1291E+02 

5 

2. 3447E+02 

2.6393E+02 

3. 4442E+02 

5. 7175E+02 

9. 4826E+02 

2. 312 IE+02 

2, S177E+02 

3. 1291E+02 

4. 2863E+02 

5.7175E+02 

2. 06G2E+02 

2.4411E+02 

2.7692E+02 

6 

5.7175E+02 

4.2863E+02 

3. 1291E+02 

2.5177E+02 

2.3121E+02 

9.4826E+02 

5. 7 175E+02 

3.4442E+02 

2. 6393E+02 

2. 3447E+02 

5. 7175E+02 

4.2883E+02 

3. 1291E+02 

7 

4.2863E+02 

5.7175E+02 

4. 2863E+02 

3. 1291E+02 

2. 6177E+02 

5. 7175E+02 

9. 4 82 6 E+ 02 

5. 7175E+02 

3. 4442E+02 

2. 6393E+02 

4. 28G3E+02 

5- 7 175E+02 

4. 2863E+02 

8 

3. 1291E+02 

4.2B63E+02 

5. 7175E+02 

4.2863E+Q2 

3. 1291E+02 

3.4442E+02 

5. 7175E+02 

9. 4826E+Q2 

5. 7175E+02 

3. 4442E+02 

3. 1291E+02 

4.2B63E+02 

5. 7 175E+02 

9 

2. 5177E+02 

3. 1291E+02 

4. 2863E+02 

5. 7 175E+02 

4.2863E+02 

2. 6393E+02 

3, 4442E+02 

5.7175E+02 

9. 4826E+02 

5. 7 175E-1-02 

2. 5177E+02 

3. 1291E+02 

4.2863E+02 

10 

2.312XE+02 

2, 5177E+02 

3. 129IE+02 

4.2863E+02 

5.7175E+02 

2. 3447E+02 

2. 6393E+02 

3.4442E+02 

5. 7175E+02 

9.4826E+02 

2.3121E402 

2. 5177E+02 

3. 1291E+02 

11 

3.4442E+02 

3. 1291E+02 

2.7692E+02 

2.4411E+02 

2.0662E+02 

5.7175E+02 

4.2863E+02 

3. 129 IE+02 

2. 5177E+02 

2.3121E+02 

9.4826E+02 

5. 7175E+02 

3. 4442E+02 

12 

3. 1291E+02 

3. 4442E+02 

3. 1291E+02 

2.7G92E+02 

2. 4411E+02 

4.2863E+02 

5.7175E+02 

4, 2863E+02 

3. 1291E+02 

2. 5177E+02 

5. 7 17 5E+02 

9. 4 826 E+ 02 

5. 7175E+02 

13 

2.7692E+02 

3. 1291E+02 

3. 4442E+02 

3. 1291E+02 

2.7692E+02 

3. 1291E+02 

4.2963E+02 

5. 7175E+02 

4. 2863E+02 

3. 1291E+02 

3. 4442E+02 

5. 7175E+02 

9.4826E+02 

14 

2.4411E+02 

2.7692E+02 

3. 1291E+Q2 

3.4442E+02 

3. J291E+02 

2. 5177E+02 

3. 1291E+02 

4.2863E+02 

5. 7175E+02 

4, 2863E+02 

2. 6393E+02 

3. 4442E+02 

5. 7175E+02 

15 

2. 0662E+02 

2. 4411E+02 

2.7692E+02 

3. 1291E+02 

3.4442E+02 

2. 3 12 IE+02 

2. 5177E+02 

3. 1291E+02 

4. 2863E+02 

5.7I75E+02 

2. 3447E+02 

2. 6393E+02 

3. 4442E+02 

16 

2. 6393E+02 

2. 5177E+02 

2.4411E+02 

2.254 1E+02 

1. 7123E+02 

3. 4442E+02 

3. 129 IE+02 

2. 7692E+02 

2. 4411E+02 

2. 0662Ei-02 

5.7175E+02 

4.2863E+02 

3. 1291E+02 

17 

2. 5177E+02 

2. 6393E+02 

2. 5177E+02 

2. 441 IE+02 

2.254 IE+02 

3. 129 IE+02 

3, 4442E+02 

3. 1291E+02 

2.7692E+02 

2.4411E+02 

4, 2863E+02 

5.7175E+D2 

4.2863E+02 

18 

2.44UE+02 

2. 5177E+02 

2, 6393E+02 

2. 5 177 E+ 02 

2.4411E+02 

2.7G92E+02 

3.1291E+02 

3.4442E+02 

3.1291E+02 

2.7692E+02 

3. 1291E+02 

4.2863E+02 

5. 7175E+02 

19 

2. 2541E+02 

2.4411E+02 

2. 5177E+02 

2.6393E+02 

2. 5177E+02 

2.4411E+02 

2.7692E+02 

3. 129 IE+02 

3.4442Ef02 

3. 1291E4-02 

2.5177E+02 

3, 1291E+02 

4.2863E+ 02 

20 

1.7123E+02 

2.2541E+02 

2.4411E+02 

2.5177E+02 

2. 6393E+02 

2. 0662E+02 

2. 441 IE+02 

2. 7692E+02 

3. 1291E+02 

3.4442E+02 

2. 3121E+02 

2 . 5177E+02 

3. 1291E+02 

21 

2.3447E+02 

2. 3 12 1E+02 

2. 0662E+02 

1.7123E+02 

2.2093E+02 

2.6393E+02 

2.5177E+02 

2. 44 1 1E+ 02 

2.2541E+02 

1. 7123E+02 

3. 4442E4-02 

3. 1291E+-02 

2. 7692E+02 

22 

2. 3121E+02 

2* 3447E+02 

2. 3 121E+02 

2. 0662E+02 

1. 7123E+02 

2. 5177E+02 

2.6393E+02 

2, 5177E+02 

2.4411E+02 

2. 2541E+02 

3. I291E+02 

3.4442E + 02 

3. 1291E+02 

23 

2. 0662E+02 

2.3121E+02 

2. 3447E+02 

2. 3121E+02 

2.0662E+02 

2.4411E+02 

2. 5177E+02 

2. 6393E+02 

2. 5177E+02 

2.44HE+02 

2.7692E+02 

3. 1291E+02 

3.4442E+02 

24 

1. 7123E+02 

2.0662E+02 

2* 3121E+02 

2. 3447E+02 

2. 3 12 IE+02 

2. 254 IE+02 

2. 441 IE+02 

2 . 5177E+02 

2. 6393E+02 

2.5177E+02 

2.4411E+02 

2.7G92E1-02 

3. 1291E+02 

25 

2.2093E+G2 

1.7123E+02 

2. 0662E+02 

2.3121E+02 

2.3447E+02 

1. 7123E+02 

2.2541E+02 

2* 441 IE+02 

2. 5177E+02 

2.6393E+02 

2.0662E+02 

2.4411E+02 

2.7692E+02 



14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

1 

2.4411E+02 

2. 0662E+02 

2. 6393E+02 

2.5177E+02 

2. 4411E+02 

2. 2541E+02 

1.7123E+02 

2. 3447E+02 

2.3121E+02 

2. 0G62E+02 

1. 7123E+02 

2. 2 093 E +02 

2 

2 . 7692E+-02 

2.4411E+02 

2. 5177E+02 

2.6393E+02 

2. 5177E+02 

2. 4411E+02 

2.254 1E+02 

2.3121E+02 

2.3447E+02 

2.3121E+02 

2. 06 62 E +02 

1. 7123E+02 

3 

3. 129 IE + 02 

2.7692E+02 

2. 44 11E+02 

2. 5177E+02 

2.6393E+02 

2. 5177E+02 

2. 4411E+02 

2.0662E+02 

2.3121E+02 

2.3447E+02 

2. 3121E+02 

2. 0662E+02 

4 

3.4442E+02 

3. 1291E+02 

2.254 1E+02 

2. 4411E+02 

2. 5177E+02 

2. 6393E+02 

2.5177E+02 

1. 7 123E+02 

2.0662E+02 

2. 3121E+02 

2. 3447 E+ 02 

2. 3121E+02 

5 

3. 1291E+02 

3. 4442E+02 

1.7133E+02 

2. 254 IE +02 

2.4411E+02 

2. 5177$+02 

2.G393E+02 

2.2093E+02 

1. 7123E+02 

2.0662E+02 

2. 3121E+02 

2. 3447E+02 

6 

2. 5177E+02 

2.3121E+02 

3.4442E+02 

3. 1291E+02 

2. 7692E+02 

2.44UE+02 

2.0GG2E+02 

2. 6393E+02 

2, 5177E+02 

2.4411E+02 

2.2541E+02 

1, 7123E+02 

7 

3. 1291E+02 

2. 5177E+02 

3. 1291E+02 

3.4442E+02 

3. 1291E+02 

2.7692E+02 

2.44 11E+02 

2. 5177E+02 

2.6393E+02 

2. 5177E+02 

2 . 4411E+02 

2.2541E+02 

8 

4. 2863E+02 

3.1291E+02 

2.7692E+02 

3. 1291E+02 

3. 4442E+02 

3. 1291E+02 

2.7692E+02 

2.4411E+02 

2. 5177E+02 

2.6393E+02 

2. 5177E+02 

2.4411E+02 

9 

5.7175E+02 

4. 2883E+02 

2.4411E+02 

2.7692E+02 

3. 1291E+02 

3. 4442E+02 

3. 1291E+02 

2.2541E+02 

2.4411E+02 

2. 5177E+02 

2. 6393E+02 

2.5177E+02 

10 

4.2863E+02 

5.7175E+02 

2. 0GG2E+O2 

2. 441 1E+02 

2.7692E+02 

3. 1291E+02 

3. 4442E+02 

1. 7123E+02 

2.2541E+02 • 

2.4411E+02 

2. 5177E+02 

2.6393E+02 

11 

2. 6393E+02 

2. 3447E+02 

5.7175E+02 

4. 2BG3E+02 

3. 1291E+02 

2.5177E+02 

2. 3121E+02 

3, 444 2 E +02 

3. 1291E+02 

2. 7692 E+ 02 

2. 44 11E+02 

2.0G62E+02 

12 

3. 4442E+02 

2.6393E+02 

4.2863E+02 

5.7175E+02 

4. 2863E+02 

3. 1291E+02 

2. 5177E+02 

3. 1291E+02 

3.4442E+02 

3. 1291E+02 

2. 7692E+02 

2.4411E+02 

13 

5. 7175E+02 

3.4442E+02 

3. 1291E+02 

4.2863E+02 

5. 7175E+02 

4.2863E+02 

3. 1291E+02 

2. 7G92E+02 

3. 1291E+02 

3.4442E+02 

3. 1291E+02 

2.7692E+02 

14 

9. 4826E+02 

5. 7175E+02 

2.5177E+02 

3. 129LE+02 

4.2863E+02 

5. 7175E+02 

4, 2863E+02 

2.4411E+02 

2. 7692E+02 

3. 1291E+02 

3.4442E+02 

3. 1291E+02 

15 

5.7175E+02 

9. 4826E+02 

2.3121E+02 

2. 5177E+02 

3. 1291E+02 

4. 28G3E+02 

5. 7175E+02 

2. 0662E+02 

2, 4411E+02 

2. 7692E+02 

3. 1291E+02 

3. 4442 E+ 02 

16 

2. 5177E+02 

2.3121E+02 

9.4826E+02 

5. 7175E+02 

3. 4442E+02 

2. 6393E+02 

2. 3447E+02 

5.7175E+02 

4. 2863E+02 

3. 1291E+02 

2. 5177E+02 

2. 3121E+02 

17 

3- 129 IE +02 

2.5177E+02 

5. 7175E+02 

9. 482GE+02 

5. 7175E+02 

3. 4442E+02 

2. 6393E+02 

4.2863E+02 

5. 7175E+02 

4.2863E+02 

3. 1291E+02 

2.5177E+02 

18 

4. 2863E+02 

3. 1291E+02 

3. 4442E+02 

5.7175E+02 

9.4826E+02 

5.7175E+02 

3. 4442E+02 

3. 1291E+02 

4. 2863E+02 

5.7175E+02 

4.2663E+02 

3, 1291E+02 

19 

5.7175E+02 

4. 2863E+02 

2.6393E+02 

3.4442E+02 

5. 7 175E+02 

9.482GE+02 

5. 7 17 5E+02 

2.5177E+02 

3. 1291E+02 

4.2863E+02 

5. 7 175E+02 

4.2863E+02 

20 

4. 2863E+02 

5.7175E+02 

2.3447E+02 

2.6393E+02 

3. 4442E+02 

5. 7175E+02 

9. 4826E+02 

2. 3121E+02 

2. 5177E+02 

3. 1291E+02 

4. 2 86 3 E +02 

5. 7 17 5E +02 

21 

2. 441 IE +02 

2. 0682E+02 

5.7175E+02 

4.2863E+02 

3.1291E+02 

2.5177E+02 

2.3121E+02 

9. 4826E+02 

5.7175E+02 

3.4442E+02 

2. 6393E+02 

2.3447E+02 

22 

2. 7692E+02 

2.4411E+02 

4.2863E+02 

5. 7175E+02 

4. 28S3E+02 

3. 1291E+02 

2. 5177E+02 

5. 7175E+02 

9.4826E+02 

5. 7 17 5E +02 

3, 4442E+Q2 

2.6393E+02 

23 

3. 1291E+02 

2.7692E+02 

3. 1291E+02 

4.2B63E+02 

5.7175E+02 

4. 2863E+02 

3. 1291E+02 

3. 4442E+02 

5.7175E+02 

9. 4826E+02 

5. 7175E+02 

3. 4442E+02 

24 

3. 4442E+02 

3. 1291E+02 

2. 5177E+02 

3. 1291E+02 

4.2863E+02 

5. 7175E+02 

4.2863E+02 

2.6393E+02 

3.4442E+02 

5. 7175E+02 

9.4826E+02 

5.7175E+02 

25 

3. 1291E+02 

3.4442E+02 

2.3121E+02 

2. 5177E+02 

3. 1291E+02 

4. 2863E+02 

5.717GE+02 

2. 3447E+02 

2. 6393E+02 

3.4442E+02 

5.7175E+02 

9. 4826E+02 



Table 4. Estimated block gravity anomalies referred to the 1967 system. 

Gravity Gravity Gravity 

1 'it. Long. Area anomaly Lat. Long. Area anomaly Lot. Long. Area anomaly 

(deg.) (deg. sq.) (mgal) No. (deg.) (deg.) (deg. eq.) (mgal) No. (deg.) (deg.) (deg. eq.) (mgal) No. 
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